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TEE INTERACTION OF BO'Jl-;DARY LAYER AND COMPRESSION SHOCK 
AND ITS EFFECT UPON AIRFOIL PRESSURE DI3TRI3UTI0JIS 

B'j H. Julian Alien, MaXe A. Heaslet, 
and G-erald Eo Nitzberg 

SUiOARY 

The mechanism of the interaction of compression shock 
Kith boundary layer is Investigated, It is shoT-Tn that the 
apparently shocliless pressure distributions observed exoeri- 
mentally at supercritical Ila.ch numbers can be accounted for 
by a marked thickening of the boundary layer for some distance 
ahead of tne shock wave. Pressure distributions with abrupt 
pressure recovery from large local supersonic Mach numbers to 
a local Mach number of about unity can be accounted for by a 
thicicenlng of the bo’undary layer over a shorter chordwise 
extent . 

At the test Reynolds number of about l,p00,000, the 
presence of an aerodynsmlcally clean or dirty surface does not 
materially influence tlie drag at high supercritical speeds. 

f 

INTRODUCTION 

DuvTing calibration of the Ames 1- by 3|-foot high-speed 
wind tunnel, a study of wa.ll interference by the image method 
was luidertalcen to determine the validity of theoretical 
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interference factors. It was found, in the covirse of the tests 
that pressure-distribution measurements at supercritical a.ir 
speeds could not alv;ays be repeated. This phenomenon was 
particularly evident at certain angles of attack. 

The fact that the phenomenon appeared only at super- 
critical speeds suggested that it was in some way brought 
about by the interaction between compression shock and the 
airfoil boundary layer. Such Interaction is possible since 
a stationary shock x^ave, v:hich only exists within a fluid 
where the velocity is above sonic, cannot penetrate to the 
sxirface of the airfoil where the veloc.lty must be zero and, 
in consequence, the pressure rise across the shock xmve must 
effect a tendency to a reversed flow xjithin the boundary 
layer . 

An Investigation x-?as undertaken to determine the cause 
of the phenomenon with a view to clarifying the present under- 
standing of compression shock and the interrelation of shock 
and boundary layer. 

SXFERIIS^TAL I^^^'ESTIG-ATION 

Ills Ames l--by 3^-f^ot high-speed wind tunnel, xvhich was 
used in these tests, is a low- turbulence wind tunnel with 
sufficient pov^er to obtain choked flow under all test condi- 
tions. The airfoil models used in the present tests \ 7 eve of 
NACA 44 i 2 airfoil section and spanned the 1-foot dimension 
of the tunnel. There 'wb.b no end leakage. One model, 
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described in reference 1, which had-a-^inch chord, was 
equipped tfltli 5 ^ pressure orifices arxd was used to determine 
the distribution of pressiure over the airfoil section. The 
other model had a 6-inch chord and was used to obtain schlieren 
photographs of the flow. 

Of the airfoil pressure distributions for the wall- 
interference study, several obtained at a Mach number of 0,75 
and an angle of attach of 4^ are of major interest for the 
present study. The Reynolds nuLiber for this test was about 
1,5 nillions. The press’Ji’e distribution for the clean airfoil 
under these conditions is shown in figure 1. 

It v;as noted during the tests that, after prolonged 
operation of the tunnel, the pressure distribution changed 
from that shown; and it was found that this change xira.s brou^t 
about by an accumulation of dirt and a pitting due to gritty 
particles striking the surface. Polishing of the airfoil 
surface permitted the original pressure distribution on the 
airfoil to be reattalned. This result suggested that the 
change in pressure distribution v;as brought about by a forward 
movement of the point of transition from laminar to turbulent 
flow. Accordingly, the pressures were measx.i.red X'^ith a thin 
spanwise strip of coarse carborundum! (No, I 6 O) fixed a.t the 
6-percent-chord station on the upper surface of the airfoil. 

An even more pronounced cliange in the pressure distribution 
of the kind previously found x-/ith the dirty surface was 
observed. This distribution is chox-Jii in figure 1. 
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The distrihution of pressures obtained vrith the clean 
airfoil surface is not coiripatible >dth the usual theory of 
supercritical floT-; since there is no sharp drop in pressure 
coefficient to indicate the presence of a shocko In order 
to study this anomaly by the schlieren raethod_, an RACA !4l|-12 
airfoil model was mounted for test on circular glass side- 
wall disks. The attachment to each disk was effected by two 
steel pins, at the 10- and lO-peroent-chord stations. In 
order to seal the end of the airfoil to the disks, thin rubber 
gaskets were used. Typical s chlieren photographs obtained 
v;ith this apparatus are shown in figure 2. It vrill be noted 
that, :juifortunately, the mounting pins and edge of the 
gaskets are visible but. they can readily be distinguished. 

In several instances, visual examination of the schlieren 
screen showed the shock wave to be blurred, indicating chord- 
wise motion of the wave. Such an unsteady flow could, of 
course, produce a pressure distribution such as that shown 
in figure 1 for the clean-surface model since the liquid 
manometer used for these tests would not respond to rapid 
changes in pressure. In many instances, however, although an 
apparently shockless pressure distribution was obtained, the 
schlieren screen indicated that a compression shock was 
present and stationary with time. This observation prompted 
an analysis of the experimental res’ults obtained. 
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ANALYTICAL IWESTIGATION 

The experimental pressure distribution obtained with a 
clean surface, shown in figure 3> is akin to the corr espondlng 
distribution of figure 1 in that, for the upper surface in 
both instances, although the llach number Is well above 
critical, evidence of shock is apparent l^r absent although the 
schlieren lehotograph of figure ^(a), corresponding to the 
pressure distribution of figure 3^ shows shocked flow. Such 
a pressure recovery, it was considered, could be produced by 
interaction of the shock \cave and boundary la^rer. 

In order to explain this, it is necessary to consider 
the general characteristics of the supersonic portions of 
pressure distributions at supercritical speeds in the absence 
of local boundary-laj^er effects. Prerdtl and Ileyer have 
shown that, when a semi-infinite uniform stream at a Mach 
namber of unity is deflected around a convex surface, the 
local Mach number attained at any point is supersonic xvith 
magnitude a function only of the total angle through vjhich 
the stream has turned. This theory is modified empirically 
in reference 2 for the fact tha.t the supersonic region of flow 
over airfoils at supercritical speeds is of only limited 
extent. It is found that the local Mach number is still 
directly related to the total angle turned through by the 
surface from the sonic point to the point under consideration. 
However, the Hach number rise for a. given angular deflection 
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is only from one-third to two-thirds that predicted by 
Prandtl and Meyer, the precise value being determined by 
the specific configuration of airfoil shape, angle of 
attack, and free-stream Mach number. 

It is knov/n from studies of boundary-layer effects at 
low speeds that the presence of a boundarj^’ layer is equiva- 
lent to a proportionate local thickening of the airfoil. 

There is no reason to believe that this effect will be 
different at high speeds. It follows that, since supersonic 
local Mach numbers, and therefore local pressure coefficients 
are directly affected by local surface contour’s, any change 
in boundary-layer thickness which markedly changes the 
effective local curvature of the airfoil surface will result 
in a marked local change in the airfoil pressure distribution 
There is an abrupt pressure increase across a shock wave 
which can be transmitted forward only through the boundary 
layer. This pressure rise prod.uces an adverse pressure 
gradient in the boundary layer which tends to thicken and, 
in some cases, actually to separate the boundary layer, 

A marked thickening of the boundary layer ahead of the shock 
wave may so modify the effective curvature of the airfoil 
surface that the original convexity, with the resultant 
falling pressure, in front of the wave can be changed to 
concavity with a resultant rising pressure in front of the 
w’ave. In other words, the effect of the thickening of the 
boundary layer can be such as to decrease the local Mach 
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numbers outside the boundary layer as the shock vrave is 
approached. The pressure rise across the shock wave is 
thereby diminished, resulting in a so-celled "softened" shock 
wave, Donaldson, in reference 3> found that the concept of the 
softened compression shock led to the successful prediction of 
the position of the terminal shock wave in supersonic nozzles. 
His prediction was based upon the assumption that this softening 
was so complete as to make the llach numbers Immediately ahead 
of and behind the terminal shock wave essentially unity. If 
the supposition is made that the pressure distribution can be 
related to the effective surface curvature and if it is ass\im.ed, 
after Donaldson, that the terminal shock wave occurs when the 
surface pressure is that corresponding to a. Hach number of 
unity, it should be possible to calcuia.te from a. Imo'wn airfoil 
pressure distribution the growth of the boundary layer e.nd the 
position of the terminal shock wave. 

In order to do this, it is necessary to relate the rise 
in pressure along the surface to the deflection in the stream 
direction. This may be accomplished by assu_ming that pressure 
changes take place on Mach lines, as was done by Ackeret and 
Busemann (reference 4, pp, 235-236), or by assuming that the 
changes occur along shock waves of Infinitesimal magnitude. 

The latter appros.ch will be adopted here. Consider nox-/ the 
uniform, two-dimensional flov; of a compressible fluid at a 
Mach number Mi > 1. If this flow is deflected through an 
angle 9 as shovi^'n in figure 5^ xvhich is less than a certain 
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prescribed vs-lue, theory indicates that a linear shock wave 
will be formed in the fluid and that the shock will pass 
through the point at which the deflection occurs. Moreover, 
classical theory has shown that it is possible to relate 
conditions on both sides of the shock wave by mes.ns of the 
equations which follow from the assumption that mass, 
momentum, and energy are conserved in ps.ssins through the 
shock (reference 4-, p. 233), As a result of these 
assumptions it can be sho^.-rn that 


P3 - = Pi V; 


3 sin 0 sin a 
cos (a-0) 


( 1 ) 


P2 - Pi = P 


V,Ysln^' a - 






1 _ '\ _2 
^ J <v+ 


7+1 


( 2 ) 


where 

p pressure 

p density 

0 angle of deflection of stream 

a angle of inclination of shock wave 

V velocity 

7 ratio of specific heats (cp/c^ = 1,4-) 

1,2 subscripts denoting conditions before and behind the 
shock, respectively. 

As 0 approaches zero, the pressure rise a.cross the 
shock approaches zero and from equation (2) it follows that 
sin a = 1 /Hi. Thus a is the Mach angle associated with 
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the free-stream Mach nximber and the shock wave becomes in the 
limit a Mach line. In the analysis that follows, an expression 
will be derived relating the change in pressure with the value 
of G . 

By setting 


^p/ci = ^ 

equations (1) and (2) may be written, I'espectivel^'', in the 
forms 


It is possible to eliminate a between the two preceding 
equations and, as a result, AP nay be expressed as a function 
of tan 0 or, inversely, the angle of deflection may be 
expressed in terns of the pressure change. Solving for tan 6 
gives 


P3 - Pi = Ap 



2 tan g tan 6 
1 + tan a tan 6 


sln*^ a ( 7 + 1 ) AP ^ 1 

)t ,r 3 
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If 0 is small, It follows that the pressure increment 
is small and, as a consequence, the right-hand member of 
equation ( 3 ) cs,n he expanded in series form. Neglecting 
powers of AP above the second yields the following 


tan 9 = AP 


2 


(AP)® 


~T~ 


(7+l)lii^ 

4(Hi^-1) 


( 4 ) 


If, in turn, AP is expressed in terms of tan 0, and 
powers of tan 0 above the second are omitted, the result is 


A.P = p Vi^ 


f 

tan 9 ^ tan^9 

( 7 + 1 ) ^ 


-1 -1 

pOif - 1 ) 




It is thus possible, under the a-ssumptions of thin— 
airfoil theory, to t reat pressure variation in a supersonic 
stream! either as a limiting case of a.ttenuated shock waves 
or by means of the previously known approach which assumes 
the pressure changes take place along Hacii lines. 

If the relation between pressure-coefficient change and 
angle-of-stream deflection given in equation [K) is used, it 
is then possible to calculate the grovrth of the boundary layer 
and the change in local llach number by a step-by-step process. 
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DISCUSSION AND CONCLUDING REiaRKS 

The- theoretical analysis was opplied to predict the 
bo’jndarv-layer growth and the formation of shock corresponding 
to some pleasured pressure distributions,, A coirparison of 
severax calculated flows with the corresponding observed flows 
is shown in figures o, 7; ^-i^d (The flows of fig, 4, as 

noted previously^ correspond to the pressure distribution of 
f ij. , 3») I^ these calcuxatxons the Lach lines corresoonding 
with the measiu’ed surface pressures have been arbitrarily 
extended linearly and are terminated at the extended Hach line 
of unity. Such a diagram clearly cannot conform with reality 
at points removed from the surface since it fails to allow for 
the proper pressure variation normal to the airfoil. In 
ps-rticular, the unity liach line, although possible close to 
the surface^ could not extend into the flow field to intersect 
the other- Mach lines as shown. In spite of these limitations, 
io is seen that, in figures 4 and 6, the terminal shock is at 
a constant chordwlse position along the span and the calculated 
and experimental flows are in good agreement, indicating that 
the previously discussed mechanism for the interaction between 
shock and boundary layer is substantially correct. The 
schlieren photographs of figures 7 and g suggest that the flow 
varies spanwise. In spite of this, reasonable agreement is 
obtained. In addition, the figures shov: that near the surface 
the actual location of the shock wave is very close to the 
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calculated position for a Mach number of unity which 
confirms Dona.ldson's hypothesis. 

With a clean surface it can be assumed that, at the 
relatively low Reynolds number of the test, the flow remained 
laminar up to the shocho Although no measurements were made 
in the nresent investigation to support this assumption, a 
similar condition vras found to exist at much higher Reynolds 
numbers in an investigation in the Ames l6-foot high-speed 
wind tunnel of a large chord NACA low— drag— oype airfoil. 

With the addition of carborundum, to the a.irfoil of the 
present investigation, tra.nsition from laminar to turbulent 
flow probably moves forwa.ra enougn so ohat a turbulent 
boundary ?.ayer exists over the region of interest. It is 
well known that while the laminar layer is easily thickened 
by an adverse pres sure gradient, the turbulent layer is more 
resistant to such a gradient. It is then to oe expected that 
the thickening of a laminar boundary layer due to the pressure 
rise across a given shock vjave will be greater tnan for a 
tiurbulent boundary layer. If it is assumed that the terminal 
shock wave occurs at a Mach number near unity, the previous 
analysis has shown that the angular deflection of the stream 
at the outside of the boundary layer has a fixed value. 

Hence the marked thickening due to the terminal shock must 
be more limited in cliordwise extent for a turbulent layer than 
for a lam.inar layer. The pressui’e distributions of figure 1 
show this to be th.e case. 
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_ 'i the cases so far considered, sufficient thickening to 
Induce shock near a Mach number of unity has been evidenced by 
the pressure distributions. It cannot be concluded, however, 
that all shocks must take place at close to unity Mach number 
since the lower-surface pressures of figure 3 would not appear 
to support such a conclusion. As a. further warning, it should 
not be concluded that whenever a pressure distribution of the 
tj^e shov.m for the upper surface in figure 3 is obtained the 
s’nock necessarily takes place near unity Mach number, unless 
the pressures are both instantaneous and persisting. Usually 
pressures are measured on licuid manometers possessing slow 
response which, in consequence, can indicate for this latter 
type of pressure distribution either a softened shock, as 
previously discussed, or a f luctuating-shock flow which in 
other instances has been observed. 

Under the. assumption that shock occurs at about the 
chordwise station vhere the local Mach number has dropped to 
unity, the pressure distributions of fig\rre 1 indicate the 
following: When the boundary- layer flov/ is changed from 

laminar to turbulent, by means of the strip of coarse carbo- 
rundum, the location of the shock v;ave is moved fon^^ard and 
the boundary— lay er separation after shock appears complete. 

As will be seen in figure 2, the growth of the boundary 
layer ahead of the terminal shockwave is usually one wherein 
the thickening- of the boundary layer is at first slow and 
increases as the terminal shock -is approached, as shown in 
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figure y{tx). The pressure distribution is then of the type 
shoT'n in figure 9(t>)» The Mach lines, which are but shock 
waves of weak strength, are generated at the boundary- layer 
s;,!rfa.ce and would form an oblioue shock of greater strength 
if they coalesced to form an envelope o However, such 
coalescence has not space enough to occur ahead of the 
terminal shock in the usual case, and hence, the oblique 
shock would not be evident. 

If, on the other hand, the boundary layer v;ere laminar 
and separated to form a X'/edge of constant angle, as shown in 
figure 10(a), the envelope of liach lines would form an 
oblique shock xi-rlxich would disappear at the terminal shock 
and the pressure distribution would be of the type shorn in 
figure lO.(b), That such types do actually occur is seen in 
figure 2 for some of the lower-surface shocks. 

It is important to note that, X'/ith the addition of the 
carborundum, the reduction of lift was sizable but that the 
relative change in drag vraa slight,, In figure 11 is shox'm 
the pressure dra,g coefficient (which at speeds x:ell in excess 
of the critical is practically the total drag coefficient) 
for the airfoil whose pressure distributions are shox\’’n in 
figure 1. It appears that, whether or not the energy loss 
in the boundary layer is increased x\>’ith a consequent reduction 
in energy loss in the shock x^^ave, the net effect on the total 
drag is practically the same. It should be noted, hox-rever, 
that in the cases considered, a softened shock is apparently 
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always attained and the conclusion majr not be iirarranted for 
other types of shocks o Finally, it should be observed that 
although the drag is relatively insensitive to changes in the 
type of boundary— layer flov; these changes produce Important 
shifts in the location of the shock wave and in the lift. 
Therefore, increase of Reynolds number, in its effect on 
thinning the boundary layer relative to the airfoil dimensions 
may have an Important effect at high liach numbers, 
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Figure 1.- Effect of surface condition on the distribution of pressures over on NACA 4412 

airfoil at 4^ angle of attack, Mach number 0.75, and Reynolds number 1.5 millions. 



Figure 3.- Distribution of pressures over an NACA 4412 airfoil section at 0° angle of attack, 
Mach number 0.772, and Reynolds number 1.5 millions. 
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Figure 2.- Representative schlieren photographs of the flow about an NACA ^12 airfoil 
at several Mach numbers and angles of attack. 
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(h) Calculated flow rait€.pn. 
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rt^uR-E 4-r Calculated and experimental flow about an nacaaaib, 

AIRFOIL SECTION AT 0°AN&LE OF ATTACK AND MACH NUMBER, 0.77^.'^ 
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Fig. 6 



(a) ExPERiMChiTAL FLOW PATTCRN. 
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AlPFOIL.SeCTfOH AT 0° ANOL€ OF ATTAC K AND MACH NUMBER. 0,J22,. 
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r«^uR.E 7»" Calculated and e>cpcp\p\eutal. flow about an naca 4Ai:l 

A! R.FOIL SECTION AT ^<*ANGl£ OF ATTACK AND MACH NUFfBFK 0.7S0> 
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Fig. 8 



Expcrimcntal FLO\x/ fattern. 
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(b) Calculated flow pattern. A-6664 
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FieuRc &.- Calculated and experimcntal flow about an nacaaaix 
AIRFOIL SECTION AT O^ANSLE OF ATTACK AND MaqH NUMBER 0,/00> 
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(i>) pressure distribution. 
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(b) PRESSURE JT/STRIBUT70N 
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FI6URE 3.- INTERACTION OF SHOCK WAVE AND BOUNDARY FIGURE 10.- INTERACTION OF SHOCK WAVE AND BOUNITARY 

LAYER FOR SINGLE-SHOCK FLOW, LAYER FOR DOUBLE-SHOCK FLOW. 
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Figure 11.- Effect of surface condition on the pressure 
drag of an NAOA 4412 airfoil at 4° angle of 

attack. 
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